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Abstract. - Starting from a microscopic model of randomly cross-linked particles with 
quenched disorder, we calculate the Laudau- Wilson free energy 5 for arbitrary cross-link densi- 
ties. Considering pure shear deformations, 5 takes the form of the elastic energy of an isotropic 
amorphous solid state, from which the shear modulus can be identified. It is found to be an 
universal quantity, not depending on any microscopic length-scales of the model. 



Introduction. - It is well known that randomly cross-linked networks of molecules un- 
dergo a transition from a fluid to an amorphous solid state, as the number of cross-links is 
increased. This sol-gel transition has been studied in detail [1-4] in recent years, and the 
structure as well as the elasticity [5, 6] is well understood. However, almost all theoretical 
studies have focussed on the near-critical region, where the analysis simplifies due to the 
existence of a small parameter - the distance to the critical point or, equivalently, the order 
parameter. In contrast, the highly cross-linked regime has hardly been investigated [7], yet it 
is particularly interesting in the application of random network models to glasses [8] . 

In this paper we consider a particularly simple network which is built from spherical 
particles, connected by harmonic springs. This allows us to access the highly cross-linked 
regime and compute the shear modulus for arbitrary cross-link concentration. We find that the 
shear modulus is independent of all microscopic length and energy scales of the model, such as 
the excluded-volume interaction, the spring constant and the length scale that characterizes 
the localisation of particles in the amorphous solid state. Instead, the shear modulus is 
completely determined by the density of cross-links or, equivalently, the gel-fraction. 

Model: randomly cross-linked particles. - We consider a system of N identical particles 
at positions {Ri,...,Rjv} in a d-dimensional volume V. Permanent cross- links connect M 
randomly chosen particles, such that a particular cross- link realization C is specified by a list 
of M pairs of monomers C = {(«i, ji), (im, Jm)}- The cross-links are modelled by harmonic 
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spnngs: 

e=l 

Here, the parameter a can be interpreted as the typical length of a cross-link. (Energies 
are measured in units of k^T.) Furthermore, a repulsive excluded-volume interaction i?ev is 
introduced to prevent the system from collapsing. The overall Hamiltonian is H = i^xunk + 
i?ev AH thermodynamic properties, including the elastic ones, can be obtained from the 
partition function Z{C) = J d'^Ri ■ ■ ■ d'^R^ which depends on the quenched disorder C 
that specifies the configuration of cross-links. 

The average over all cross-link configurations C is called the disorder average and denoted 
by [•]. A realistic choice of the cross-link distribution is the Deam- Edwards distribution [1]: 

1 fV 



Its parameter /x^ determines the probability for a cross-link to be formed, and therefore con- 
trols the sol-gel transition. The average cross-link density is given by [M]/N = and 
the standard deviation, relative to the mean, vanishes in the thermodynamic limit. Thus, 
/j,^ = 2[M]/N is the average coordination number, i.e., the average number of particles to 
which a certain particle is connected. For a detailed discussion of the Deam-Edwards distri- 
bution, see [3,7]. 

Order parameter and free energy. With help of the replica technique we obtain the 
disorder-averaged free energy = — [InZ] = — lim„^o(-Zn-i-i — -2i)/(n-Ei). The replica par- 
tition function 2n+i = J 'Dfle~^" is represented as a functional integral over the order- 
parameter field, whose expectation value 

i=l 

quantifies the probability of finding a particle at position x*^*^) in replica 0, at x*^^^ in replica 
1, etc. Here, the average {■)sn is evaluated with the statistical weight e""^'^. To simphfy 
the notation, we have introduced d(n -I- l)-dimensional hatted vectors x = (x(°\ x'")). 
Previous work [9] has shown that the saddle-point solution for the order parameter has the 
following simple form: A fraction of monomers (1 — Q) is delocalized and can be found 
anywhere in the sample with equal probability. The remaining fraction Q is the "infinite" 
cluster and each particle i performs Gaussian fluctuations in each replica about a mean 
position Yi with localization length ^j. The mean positions are randomly distributed and 
the localization lengths ^ follow the distribution V{^^). The resulting order parameter is 
explicitly given by 



^(xa-y)^ 



2^2 



(4) 



with self-consistent equations for Q and 'P(^^) given in [7,9]. Note that the order parameter 
does not depend on YlZ=o 

The Landau- Wilson free energy 

,,2 2 / f r\d{n+l) r> / 
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is particularly simple, because we consider cross-linked particles instead of cross-linked chains. 
This will allow us to perform computations in the highly cross-linked limit. The function 
A(A;) = exp(—ia^A:^) reflects the harmonic potential for the cross-links Q), and fc :— (k'"' , k^")). 
The Symbol means summation over all k with at least two non-zero components k^") . 

Spontaneously broken translational invariance. - The effective replicated Hamiltonian 
is invariant under uniform translations of all particles, separately in each replica, i.e., the 
transformation Rf Hf + u" is an exact symmetry of the Hamiltonian(^). In the amorphous 
solid phase, which is characterized by random localisation of a finite fraction of particles, 
this symmetry is spontaneously broken, and the only symmetries remaining are common 
displacements of all replicas, i.e., Rf Rf + u, reflecting the macroscopic translational 
invariance of the amorphous solid phase. (For a detailed discussion see [5,6].) Consequently 
there is a whole family of order parameters that each give rise to the same free energy and 
are related by replica-dependent, spatially uniform translations. To distinguish common and 
relative displacements of the replicas, it is helpful to split up the d{n + l)-dimensional hatted 
vectors x into d-dimensional longitudinal ones xy = ^^^^ Sq=o -^^""^ ^^'^ dn-dimensional 
transverse ones x± = x — ^n+i ^"^ll ' ■•■'■X||)- The manifold of symmetry-related saddle points 
is parameterized by mj^ = ii — --^===(u|| , U|| ): 

The situation is illustrated in fig. where the order parameter is shown for two replicas in 
one space-dimension. The order parameter does not depend on X|| , so that it is represented 
by a rectilinear ridge having its maximum height along the axis x± = and having width 
^. The family of symmetry-related solutions is generated by rigidly translating the ridge to a 
new position: xj_ = uj_. (The delocalized particles contribute only an unimportant constant 
background (1 — Q)/V"^^, which has been ignored in the figure.) 

Goldstone fluctuations. - Low energy excitations are constructed by particle displace- 
ments u_l(x||) that depend on the position in the sample [5,6]. The ridge is not translated 
rigidly, but deformed by the position-dependent vector m_l(x|[), as indicated in fig. These 
excitations do cost energy; only in the long wavelength limit do we expect the energy of ex- 
citation to vanish with an inverse power of the wavelength. The generalized order parameter 
fluctuations take the form [10]: 

We only consider pure shear deformations, and hence require that |det((5iy_p -f- d^u'^)\ = 1. 
This constraint guarantees the preservation of volume, and simplifies for small deformations 
to di,u'^ = 0. 

We insert the Ansatz {TJ into the general expression (0 for the free energy and only con- 
sider small distortions, so that it is sufhcient to keep uj_ only to the lowest order, (9u_l/9x|| )2. 
Higher orders in u±^ and higher-order derivatives are neglected. Restoring units of energy, the 

{}) The effective replicated Hamiltonian is also invariant under uniform rotations, separately in each replica. 
However the spontaneous breaking of this symmetry does not give rise to Goldstone modes and is not discussed 
further here. 
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Fig. 1 - Ansatz for the order parameter f2(a;'^', a;'^') in real space for two replicas in one dimension. 
White indicates high values of Q. To simplify, we have assumed a fixed localization length in this 
figure. Q{x^^\x^^^) is the probability of finding a particle at position in replica 1 and at in 
replica 2. For the rotated coordinate system a;|| = {x^^^ + x^^^)/\/2 and x± = {x^^'' — x^^^)/2. In (a) 
the displacement vector u± = 0, and in (b) u± = sina;||. (Note, however, that for physical systems 
typically ^ <^ wavelength of itx.) 



resulting free energy is 

with the shear modulus G given by 

G = (^ii^Q - 1 + exp(-/z2Q) - no k^T , (8b) 

where no = N/V is the mean particle-density. 

The part of the Landau- Wilson free energy H8a|l that does not depend on u± determines the 
gel fraction and the distribution of localization lengths that was discussed in [7]. In particular, 
differentiating Su^=q with respect to Q yields a simple relation between the gel fraction Q 
and the average coordination number /x^: exp(— /x^Q) = 1 — Q, which is shown in fig. EJ It 
was previously derived in many other contexts [3,7,9]. 

Shear modulus. - In eq. (j8b|) . we identified the shear modulus G in the free energy. As 
one would expect, it vanishes in the liquid phase {p^ < !)■ In the solid phase {p? > 1), it 
increases linearly with particle density no and depends only on the cross-link density, given 
the gel fraction Q = Q{^j,^). 

Close to the gel transition /i^ > 1, the shear modulus can be expanded in terms of e := /i^ — 1, 
resulting in power-law behavior (see fig. |3Ji) : 



G = ^nofcBre3 + 0(e4) 



(9) 
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Fig. 2 - Fraction of localized particles Q, depending on the average coordination number. The infinite 
cluster forms at the critical value /i^ = 1- 



For the highly-connected regime, /u^ 1, G reduces to (see fig.|3|D): 



(10) 



The scaling G ^ close to the critical point has been derived previously [5,6]; the 
critical exponent is a result of mean field theory and expected to be modified by interacting 
fluctuations below d = 6 [11]. In the highly cross-linked limit, fluctuations are presumably 
weak so that our result is expected to hold even beyond the Gaussian expansion around 
mean field theory. The result we get is reminiscent of the classical theory of rubber elasticity, 
G = Ugk^T , with the density of strands rig being replaced by the density of cross-links: 
y'^no = [M]/V. 

In order to compare the calculated shear modulus with experimental results, note that, if 
the cross-link length is small compared to the particle size, it is unlikely for one monomer in 
d = 3 (d = 2) dimensions to be connected to more than 12 (6) other monomers. Therefore, 
average coordination numbers jj? > 12 (> 6) are inaccessible in these systems. Only if long 
tethers can form, reaching further than the nearest neighbour, higher average coordination 
numbers /i^ are possible. 

Conclusions. - We have determined the shear modulus G of a randomly cross-linked sys- 
tem. It only depends on the particle density no = N /V and the average coordination number 
/i^. Interestingly, neither the distribution of localization lengths 'P(C^), which entered the or- 
der parameter in eq. Q, nor the typical cross-link length a influence the macroscopic behavior 
of the material. This result is consistent with the observation that the shear modulus deter- 
mines the response of the system to a shear deformation on the longest scales: The wavelength 
of the shear deformation has to be larger than all microscopic length scales. Furthermore, one 
observes that G ~ T. Therefore, the shear elasticity is purely entropic. In contrast, the bulk 
modulus should depend on the strength of the excluded-volume interaction, and hence it is 
not expected to be as universal as the shear modulus. 

We also stress that the shear modulus does not depend on the spatial dimension d, apart 
from the fact that the maximal coordination numbers depend on dimensionality and determine 
the physical range of . Again, this could be expected, because we only consider a Gaussian 
expansion around the saddle point; interactions between fluctuations will most likely cause a 
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Fig. 3 - Shear modulus G divided by the particle density no and k^T , depending on the average 
coordination number ^ . The solid line shows the exact expression ISbll . the dashed lines are the 
approximations valid (a) close to the sol-gel transition, eq. lO, and (b) for highly connected systems, 
eq. GOJ. 

dependence on dimensionality. 

It is also of interest to study the local elasticity, i.e., the response to shear deformations of 
finite or even small wavelength. Random networks are strongly inhomogeneous, so that one 
expects the elasticity to fluctuate spatially. Work along these lines is in progress. 
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